Interplay between symmetry and spin-orbit coupling in graphene nanoribbons 
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We study the electronic structure of chiral and achiral graphene nanoribbons with symmetric 
edges, including curvature and spin-orbit effects. Curved ribbons show spin-split bands, whereas 
flat ribbons present spin-degenerate bands. We show that this effect is due to the breaking of spatial 
inversion symmetry in curved graphene nanoribbons, while ffat ribbons with symmetric edges possess 
an inversion center, regardless of their having chiral or achiral edges. We find an enhanced edge-edge 
coupling and a substantial gap in narrow chiral nanoribbons, which is not present in zigzag ribbons 
of similar width. We attribute these size effects to the mixing of the sublattices imposed by the 
edge geometry, yielding a distinct behavior of chiral ribbons from those with pure zigzag edges. 
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I. INTRODUCTION 

The crucial interplay between structure and electronic 
properties of graphene is among the most attractive fea- 
tures of its derived nanomaterials. Both carbon nan- 
otubes (CNTs) and graphene nanoribbons (GNRs) show 
promising characteristics for spintronic devices. Recent 
progress in experimental techniques has allowed for the 
fabrication of graphene nanostripes by using electron- 
beam lithograph}/^ or by unrolling CNTs.^ These ribbons 
could be used in electronic devices, such as field-effect 
transistors (FETs),^ opening new perspectives for nano- 
electronics. 

The presence of localized edge-states in GNRs, theo- 
retically predicted,'^^ and experimentally proved,^ con- 
fers them distinct properties. GNRs have attracted a 
great amount of theoretical work, but mostly focused 
on high- symmetry zigzag and armchair achiral ribbons. 
Zigzag ribbons have edges states with different spin po- 
larizations, while armchair nanoribbons do not have edge 
states. The edges of minimal^ chiral ribbons can be con- 
sidered as a mixture of armchair and zigzag edges, thus 
having edge-localized states stemming from their zigzag 
part. Although the evolution of the nanoribbon band 
structure upon the change of chirality has been recently 
addressed,^ "^^ these systems have been nonetheless much 
less studied and many aspects remain to be explored. 

The seminal work of Kane and Mele^ triggered the 
interest on new quantum phases of matter and on the 
spin-orbit coupling (SOC) effect which, although known 
to be small in graphene, gives rise to important physics. 
In particular, the quantum spin Hall (QSH) phase has 
been widely addressed. ^I^nil 

Curvature is known to enhance spin-orbit interaction; 
its importance in SOC effects has been theoretically 
investigated for the honeycomb lattice, especially for 
Q]\jrpg[T5H20| ^^^ experimentally confirmed. ^^ Hybridiza- 
tion of TT and a orbitals, decoupled in fiat graphene. 



is enhanced by curvature and thus SOC effects are 
bolstered. ^^ The interplay between curvature and SOC 
in GNRs has been mostly focused in achiral ribbons, 
with highly symmetric zigzag and armchair edges. ^^^^ 
For zigzag GNRs, dispersionless edge bands in the fiat ge- 
ometry were found to become dispersive because of SOC 
effects. Both n and a edge states remain spin-filtered 
in the curved geometry, still localized at the boundaries 
of the ribbon, albeit with an in-plane spin component 
and a localization length larger than for the fiat case.^ 
Recent experiment j^ * ^^ * ^ ^ on chiral GNRs obtained by 
unzipping CNTs show a reminiscent curvature. Scan- 
ning tunneling spectroscopy measurements revealed the 
presence of one-dimensional edge states, with an energy 
splitting dependent on the width of the ribbon. ^^ By com- 
parison with calculations employing a 7r-band model with 
a Hubbard term the width dependence of the edge state 
gap was in terpre ted as a consequence of spin-polarized 
edge states .1^^^ Hence, the study of curvature effects in 
GNRs is relevant from the experimental and theoretical 
viewpoint. 

In this paper we address the study of this ampler class 
of ribbons with chiral edges, focusing on the differences 
of SOC effects in fiat and curved nanoribbons. We sum- 
marize our main results. 

(i) We find that the bands of both, chiral and achiral, fiat 
ribbons with symmetric edges are at least two-fold spin- 
degenerate due to spatial inversion symmetry. Curving 
the ribbons breaks this symmetry, thus yielding spin-split 
bands except for the time-reversal protected special sym- 
metry k points. 

(ii) We find a gap in all chiral ribbons, despite they have 
a zigzag edge component. Boundary conditions in chi- 
ral ribbons mix both sublattices at each edge. This en- 
hances edge-edge coupling, which results in a substantial 
gap without invoking electron interactions, 
(iii) Curvature augments spin-orbit efiFects in GNRs, 
yielding a larger splitting in the spin-split bands. In fact. 



curvature may induce metallicity in ribbons which have 
a gap in the planar form. 

(iv) The spatial distribution of edge states depends on 
curvature and chirality. While zigzag ribbons are known 
to have spin-filtered states at the edges, in narrow {^ 40 
A) chiral ribbons edge states can have nonzero density 
at both edges, due to the edge-edge coupling. This size 
effect is more evident in ribbons with chiral angle close 
to 30°, i.e., that of the armchair edge, for which the sub- 
lattice mixing is stronger. 

This paper is outlined as follows. Section [ll| describes 
the structure and geometry of the ribbons studied. Sec- 



T=Tz+2Ta 



tion III gives some symmetry considerations concerning 
the role of spatial inversion in flat and curved general rib- 
bons. Section IV contains the description of the model 
Hamiltonian and calculation method. Section V presents 
the results, including spin orbit interaction, for ribbons 
of different width in flat and curved geometries. Finally, 
in Section VI we discuss our results and final conclusions 
are drawn. 



II. GEOMETRY 

We focus on chiral ribbons with symmetric minimal 
edges ,1^ obtained from unrolling chiral carbon nanotubes. 
The ribbon is thus characterized by the edge vector 
T = nai + ma2, where ai and a2 are the primitive 
vectors of graphene, and the width vector W. The 
widths considered are therefore given by an integer mul- 
tiple of H^ defined as the smallest graphene lattice vec- 
tor perpendicular to T, as depicted in Fig{l] For a 
given T, H is uniquely determined up to a global ±1 
factor. As W = MH^ we will denote the ribbons by 
M{n,m)^ where M states the width of the ribbon, and 
(n, m) indicates the minimal edge. All minimal edges 
can be decomposed in a zigzag and an armchair part,!^ 
T = nzTz + uaTa, with Tz = ai and T^ = ai + as- 

The chirality of the ribbon is specified by the chiral an- 
gle between the translation vector (n, m) which defines 
the edges and the zigzag direction (1,0). 

We take into account different curvatures in the 
transversal direction for a given flat GNR, with no 
stretching allowed. Curvature is denoted by the angle cp 
spanned by the ribbon from its curvature center, ranging 
from for a flat ribbon to a value of 27r, which corre- 
sponds to a nanotube with cut bonds along its length. 
The degree of curvature is controlled by the angle and by 
the diameter of the cylindrical configuration. 



III. SYMMETRY CONSIDERATIONS 

Carbon nanotubes are classified as achiral and chi- 
ral according to their having a symmorphic or non- 
symmorphic symmetry group respectively. This means 
that chiral tubes possess an spiral symmetry, so that 
there are two enantiomers for each chirality, while achiral 




FIG. 1. (Color online) Geometry of the 2(3,2) GNR high- 
lighted in dark gray on a graphene sheet, showing its trans- 
lation vector T = Tz + 2Ta and its width vector W = 2H, 
where H = — 7ai + la2. The unit cells spanned by T and H 
or W are indicated with dotted lines. 



tubes are equal to their mirror image, i.e., achiral tubes 
present space inversion symmetry while chiral tubes do 
not.^^ ^^ Graphene nanoribbons, as their siblings carbon 
nanotubes, are customarily classified as achiral and chi- 
ral according to their edge shapes. In GNRs this clas- 
sification is related to the crystallographic orientation of 
the boundaries: ribbons with zigzag and armchair edges 
(derived from armchair and zigzag CNTs, respectively) 
are called achiral, and those obtained from chiral tubes 
are called chiral GNRs. However, these so-called chiral 
ribbons with symmetric edges do have an inversion cen- 
ter. Upon bending the ribbon the inversion symmetry 
is lost. This feature is crucial when considering SOC ef- 
fects. Fig. |2] shows an example of a flat (left panel) and 
curved (right panel) unit cell of the (3,2) ribbon, the lat- 
ter with cp = 7T. A symmetry center is indicated in the 
planar geometry. No such an inversion center exists in 



the curved nanoribbon. 




FIG. 2. (Color online) Schematic geometry of a flat (left) and 
curved (right) 1(3,2) GNR. The curvature angle is (f = tt. 
Two pairs of equivalent atoms under spatial inversion sym- 
metry are highlighted, and the inversion center marked as O. 



portant contribution of the crystal potential to the spin- 
orbit coupling is close to the cores, the Hso contribution 
takes the form 



Hso = Y1 



h 1 dVi 
Am?c^ Vi dvi 



L-S = XL-S, 



where the spherical symmetry of the atomic potential Vi 
has been assumed and r^ is the radial coordinate with 
origin at the i atom. L stands for the orbital angular 
momentum of the electron, and S is the spin operator. 
The parameter A is a renormalized atomic SOC constant 
which depends on the orbital angular momentum. Notice 
that the Hso terms only couple p orbitals in the same 
atom. Considering the spin parts of the wave functions, 
the Hamiltonian matrix has SNa x SNa elements, Na be- 
ing the number of the C atoms in the unit cell of the 
GNR and 8 corresponds to the four orbitals per spin of 
the sp^ basis set. The total Hamiltonian in the 2x2 
block spinor structure is given by 



Notice that a flat nanoribbon with different edges lacks 
inversion symmetry. This situation is of experimen- 
tal interest: most likely, actual ribbons will not have 
symmetric edges. Such asymmetry can be achieved ei- 
ther by adding or removing atoms to an originally edge- 
symmetric ribbon, or by altering the bond lengths in one 
of the edges. Dissimilar bond lengths may arise as a re- 
sult of a different functionalization on the two edges of 
the nanoribbon. 



IV. THEORETICAL MODEL AND 
COMPUTATIONAL METHODS 

We calculate the band structure of graphene within the 
empirical tight-binding (ETB) approximation. Although 
the TT-orbital tight-binding model is known to capture the 
low energy physics of graphene, since we are interested in 
SOC effects, we consider here an orthogonal four-orbital 
25, 2px^ 2py^ 2pz basis set. This allows for the inclu- 
sion of the intrinsic SO terms within the conventional 
on-site approach. The matrix Hamiltonian is built fol- 
lowing the Slater-Koster (SK) formalism up to nearest- 
neighbor hopping. We use the parametrization obtained 
by Tomanek-Louie for graphite. ^^ The expression of the 
one-electron ETB Hamiltonian is 



^o = E^-+ E 






H.C., 



(1) 



(i,J>,/3,^ 



where Ca represents the atomic energy of the orbital a, 
(z, j) stands for all the atomic sites of the unit cell of the 
GNR, and c^j" and cf^ are the creation and annihilation 
operators of one electron at site z, orbital a and spin 
5, respectively. We focus on neutral graphene, thus no 
doping effects are addressed. 

SOC effects are included by adding an atomic-like term 
Hso to the Hq Hamiltonian. Assuming that the most im- 
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The total Hamiltonian H incorporates both, spin con- 
servation and spin- flip terms. The spin-conserving diago- 
nal terms act as an effective Zeeman field producing gaps 
at the K and K' points of the graphene BZ, with oppo- 
site signs .^ By exact diagonalization of the matrix H 
we obtain the band structure of GNRs. As explained in 
the previous Section, the curved geometry is obtained by 
isotropically bending the ribbon in the width direction, 
without changing the distance along its length. Thus, no 
bond stretching is included along the ribbon axis. We 
do not consider reconstruction or relaxation of the edges, 
nor passivation of the dangling bonds. 

The value of the SOC constant for C-based materials 
is not well established and it is still under debate. Some 
theoretical estimates gave A = 1/ieV for graphene ,1221211 
much smaller than the atomic SO coupling, 8 meV. Tak- 
ing into account the role of d orbitals, this value raises 
to A = 25/ieV.^^ Accurate measurements of SOC are 
difficult to perform in graphene because external effects 
such as substrates, electric fields, or impurities may mask 
its value. However, recent experiments in CNT quan- 
tum dots have reported spin-orbit splittings substantially 
higher than those theoretically predicted: Kuemmeth et 
a/.^give a maximum splitting of 0.37 meV in a CNT of 
diameter 5 nm. One possible explanation for this energy 
splitting is that a higher value of A should be consid- 
ered; as indicated by Izumida et al^^^ those measure- 
ments are compatible with A = 14 meV. More recently, 
Steele et al.^^ have presented evidence of large spin-orbit 
coupling in CNTs, up to 3.4 meV, an order of magnitude 
larger than previously measured and the largest theoret- 
ical estimates. Furthermore, transport experiments re- 
port spin-relaxation times in graphene 1000 times lower 
than predicted. ^^ This is compatible with a larger value of 



the SOC coupling than those given by previous theoret- 
ical estimates. ^^^^ Although small, its effects in GNRs 
could have important consequences when considering the 
spin degree of freedom, as it has been experimentally 
shown in CNTs.^^ ^^ For the sake of clarity, we choose 
for the Figures a spin-orbit interaction parameter A = 0.2 
eV. 

The spin-orbit contribution to the Hamiltonian, Hso 
is linear on A. We have checked that, for small values 
of this parameter, such as the one employed here and 
those of physical relevance, the eigenvalues of the full 
Hamiltonian H are basically a linear correction to those 
without the SOC term, Hq. Therefore, the spin-orbit 
splittings are proportional to A and the results presented 
in this work can be scaled accordingly. 



V. RESULTS 

We have calculated the electronic properties of many 
different chiral GNRs, varying their width and curva- 
ture. All calculations have been performed with the four- 
orbital parameterization explained in Sec. IV There- 
fore, they show some differences with respect to the 
widely- used one-orbital approach."^ We present here the 



band structures for three representative chiral ribbons 
(M(7, 1), M(5,2), and M(3,2)) and for comparison, 
some zigzag GNRs of different widths. The M(7, 1) rib- 
bon has a chiral angle = 6.58°, close to the zigzag 
direction; the M(5,2) has 6 = 16.02°, and the M(3, 2) 
GNR has 6 = 23.41°, closer to the armchair direction. 
The M(7, 1) and M(3, 2) GNRs have the same unit ceh 
with 76 atoms, but with different orientation, i.e., the 
H and T vectors are interchanged. The M(5, 2) ribbon, 
in the intermediate chirality range, has a unit cell with 
52 atoms. The (5,2) edge has three armchair {A) and 
two zigzag (Z) units, so for the infinite system there are 
two possible arrangements of the armchair and zigzag 
units with the same edge vector. We choose the one with 
all zigzag units together, the ZZZAA. While the edge 
states in a semi-infinite graphene sheet or in very wide 
ribbons are the same irrespectively of the sequence, for 
narrow ribbons some differences in the band structure of 
distinct edge arrangements may arise. 



A. SOC and inversion symmetry 

1. Curvature effects 

In systems with time reversal and spatial inversion 
symmetry, spin-orbit interaction does not lift spin de- 
generacy, according to Kramers' theorem. If spatial in- 
version symmetry is not present, the states are spin- split 
except in the k points protected by time reversal invari- 
ance. In symmetric-edge GNRs, spatial inversion sym- 
metry is broken by curving the ribbon, as indicated in 



Fig. |2j The importance of the broken inversion symme- 
try is shown in Fig. |3) where the electronic structures 
of the 1(7,1) ribbon calculated with and without SOC 
terms are depicted for (a) the planar configuration and 
(b) the curved one with (^ = 27r. This angle corresponds 
to a maximally curved geometry without overlapping the 
edges of the ribbon; it is equivalent to an open carbon 
nanotube with circumference equal to the width of the 
GNR. SOC effects are clear: all degeneracies, including 
spin, are lifted in the curved ribbon (b), while in the fiat 
system (a) the bands remain spin-degenerated. The only 
noticeable difference in the flat case is a small shift in 
the bands at the F point, as it can be observed in the 
zoom of Fig. |3] Otherwise, the effect of SOC is negli- 
gible. However, a large splitting is found in the curved 
ribbons, greater for the conduction bands. This is due to 
the interaction of edge states with higher-energy bands 
(see zoom in Fig. |3] (b)), which in fact are completely 
hybridized due to curvature. Throughout most of the 
Brillouin zone (BZ), one of the spin-split bands has an 
upward shift in energy, whereas the other band under- 
goes a downshift. Thus, the bands without SOC mostly 
lie between the spin-orbit-split bands, as it can be seen 
if Fig. [3](b), especially in the zoom. 

Notice that GNRs can be made metallic because of 
the curvature, as seen in Fig. [3J the gap observed in 
the flat 1(7,1) ribbon (panel (a)) is still present at F in 
the curved (cp = 27r) geometry, but this latter case the 
ribbon is metallic due to the band bending produced by 
curvature-induced hybridization (panel (b)). 

Although the spin is no longer a good quantum num- 
ber, the expectation value of the spin operator shows 
that in the flat geometry the total spin is normal to the 
ribbon. Curvature provokes the appearance of a small 
component in the in-plane direction which increases with 
curvature. 

Fig. |4] illustrates the effect of curvature. It presents the 
band structure of the 1(3,2) ribbon for three bending an- 
gles, namely, (a) 0.57r; (b) 1.27r, and (c) I.Stt, with (black 
dots) and without (red/gray dots) SOC. The curved ge- 
ometries are shown above each band panel. In a wide 
M(3, 2) ribbon there are four edge bands at eV ex- 
tending from |rX to X.^ In the case depicted in Fig. 
|4] a large gap opens between the occupied and unoccu- 
pied edge bands due to size effects, which we discuss later 
on. There is a general increase of the band splitting with 
growing curvature, as expected, due to the increment of 
the cr-TT hybridization produce d in t he curved ribbons, 
analogous to the effect predictecP^EU and experimentally 
measured in CNTs^. Moreover, band splitting in GNRs 
is anisotropic, band- and /c-dependent, as also found in 
CNTs. 



2. Edge modification 

As discussed above, planar GNRs with dissimilar edges 
also lack inversion symmetry. Different edges can be 




FIG. 3. (Color online) Electronic structure of the 1(7,1) GNR 
calculated with (black dots) and without (red/gray dots) 
SOC, considering (a) the four-orbital sp^ basis set in a planar 
geometry; and (b) the same sp^ basis but for a curved geome- 
try with (/? = 27r. Zooms of the edge bands are included close 
to each panel. 




FIG. 4. (Color onhne) Electronic structure of the 1(3,2) GNR 
with different curvatures: (a) ip = O.Stt (b) (/? = tt, and (c) 
(p — I.Stt. The bands with SOC are shown in black; bands 
without SOC are in red (gray). 



achieved either by adding or removing atoms in a sym- 
metric edge ribbon, or by altering the bond lengths at 
one edge of a symmetric GNR. We have explored the 
magnitude of this effect, calculating the changes in the 
band structure of a planar 1(5,2) modified ribbon. We 
have considered two types of modifications: two atoms 
of one edge have been removed, and the bond length 
of the zigzag atoms at one edge have been changed 10%. 
Both modifications break the spin degeneracy of the band 
structure; however, the splitting is 2 orders of magnitude 
smaller than the that achieved by the effect of curvature. 



Thus, in what follows, we concentrate on the curvature 
mechanism as a means to break inversion symmetry. 



B. SOC effects and Avidth of the ribbon 

In order to explore the interplay of SOC and the width 
of the ribbon we consider first flat GNRs, since in curved 
geometries they could be masked by other effects, such 
as hybridization. We have performed calculations for dif- 
ferent chiralities, verifying that there is a gap in planar 
chiral nanoribbons that decreases with increasing width. 
Panels (a) to (c) of Fig. [s] demonstrate this effect for the 
M(5,2) GNRs. The 1(5,2) ribbon has a substantial gap, 
around 0.4 eV, while the flat bands around Ep for the 
4(5,2) ribbon in panel (c) are clearly identified as edge 
bands for their dispersionless character near the Brillouin 
zone (BZ) boundary X. Nonetheless, the gap can be dis- 
cerned in the inset of Fig. p](c), as stated above. 

Comparison with high-symmetry zigzag ribbons of 
similar widths shows a striking difference. Figure [5] 
demonstrates that zigzag ribbons have a negligible en- 
ergy gap, around 0.1 meV for the narrowest case depicted 
(Fig. [H] (d), W = 8.52 A), while for the 1(5,2) ribbon of 
similar width (8.87 A) the gap is around 0.4 eV. 

The gaps between edge bands in chiral ribbons are due 
to the stronger coupling between edge states localized at 
each boundary. In zigzag ribbons the atoms at oppo- 
site edges belong to different sublattices, whilst in chiral 
ribbons boundary conditions at each edge mix the two 
sublattices, coupling the states located at the two edges. 
This results in a band gap due to quantum size effects, 
without invoking electron-electron interactions. 

In flat ribbons, SOC induces a tiny shift of the energy 
bands; in zigzag ribbons, it turns the flat edge bands into 
dispersive ones. As inversion symmetry is preserved, all 
bands remain spin degenerated. 

Now, we include curvature in order to enhance SOC 
effects and break spin degeneracy, as discussed in the 
previous Section. For narrow chiral ribbons the gap due 
to quantum size effects is rather large, as it can be seen 
in Figs. [5] (a) and (b), so the effect of SOC is the afore- 
mentioned energy shift and, most importantly, the spin 
splitting of the bands. We focus on the widest ribbons, 
namely, the 4(5,2) and the 8(1,0), with a smaller quan- 
tum size gap, and consider the same curvature radius for 
both ribbons, R = 6.274 A, which yields a angle cp = 1.87r 
and 1.737r respectively. For these cases, the bands closer 
to Ep are strongly deformed. These happen to be edge 
states, so its behavior gives rise to a more interesting sit- 
uation than in the large gap ribbons, as it is illustrated 
in Fig.jGJ 

Fig. [6] (a) shows the zigzag case, with a noticeable dis- 
persion in the edge bands. The zoom shows that the 
bands with SOC are spin-split, with a crossing point 
slightly shifted with respect to that of the bands with- 
out SOC. Fig. [6] (b) shows the chiral 4(5,2) GNR; 
here besides the energy shift and spin-splitting of the 
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FIG. 5. (Color online) Band structures for flat ribbons cal- 
culated with the sp^ basis set, including the SOC term for 
the following GNRs: for (a) 1(5,2); (b) 2(5,2); (c) 4(5,2); (d) 
2(1,0); (e) 4(1,0); and (e) 8(1,0). Insets show zooms of the 
bands near Ef- 



SOC bands, there is a slight displacement of the Fermi 
wave vector, which is no longer at X Although small, 
SOC effects have important consequences for the trans- 
port properties of curved GNRs: spin- filtered channels 
arise due to the interplay of SOC and curvature, and 
for wider ribbons, even chiral GNRs present these spin- 
filtered channels in the low-energy region. 



C. Chirality and spatial distribution of edge states 

Edge states are among the most important features of 
GNRs.^ ^^ It is interesting to explore how chirality affects 
the behavior of these states. For the sake of simplicity, 
we focus on flat geometries; generalization to curved ge- 
ometries is straightforward. In Fig. [7] we present the 
electronic densities of edge states belonging to flat rib- 
bons with different chiralities. The two states chosen 
correspond to the edge bands closer to the Fermi level, 
near the high-symmetry point to which the edge bands 
converge for large widths. Opposite spin polarizations 
are indicated with distinct colors; since we are dealing 
with fiat ribbons, the spin direction is perpendicular to 
the plane of the ribbon. We choose a k value slightly 
displaced from the symmetry point (either F or the BZ 
boundary X) in order to avoid degeneracy due to time re- 
versal invariance, and ribbons of similar width in order to 




FIG. 6. (Color online) Band structures for curved ribbons 
calculated with the sp^ basis set, for the following GNRs: (a) 
8(1,0); (b) 4(5,2). The bands with SOC are shown in black; 
bands without SOC are in red (gray). Insets show zooms of 
the bands near Ef- 



compare size effects. In pure zigzag ribbons, edge states 
are spin-filtered: each edge state has a well-defined spin 
orientation and it is located at one edge. This is illus- 
trated in Fig. [7|(a), which shows the square modulus of 
the wave functions for two edge states corresponding to 
a 10(1,0) zigzag nanoribbon, of width equal to 42.61 A. 
These two edge states of zigzag ribbons live in opposite 
sublattices, and their probability density is mostly con- 
fined to the atoms with coordination number two which 
constitute the geometrical edge. Fig. [t] (b) shows the 
density for two edge states of a 4(7,1) ribbon of width 
42.89 A, close in chirality to the zigzag case. Similarly 
to the zigzag ribbon, the density of each edge is mostly 
located in one sublattice. It is not homogeneously dis- 
tributed, being mostly in atoms with coordination num- 
ber two, although there is some appreciable weight in 
inner atoms close to the edge. For the 1(7,2) GNR (of 
width 34.88 A and i9 = 12.21°), Fig. [t] (c), which has 
a larger chirality angle, the wave function extends more 
into inner atoms, especially close to the armchair part of 
the boundary. 

For chiralities closer to the armchair, the states have a 
nonzero density at both edges: panel (d) shows the two 
states close to Ep for a 4(5,2) ribbon of width 35.48 A. 
These two states live at both edges simultaneously, with 
more inner atoms with nonzero density. This is more dra- 
matic in panel (e) of Fig. I7J which depicts the edge states 
for a 2(3,2) ribbon of width 37.40 A. As its chirality is 
closer to the armchair case, edge states have a greater 
penetration into the inner part of the ribbon. In order 
to localize these states and obtain the quantum spin Hall 
phase, a larger value of the SOC constant (A ~ 4 eV) 
is needed. This behavior has been found for armchair 
ribbons within the Dirac model^^ and the ETB model.^ 
The SOC strength required increases with the chiral an- 
gle and decreases with the width of the ribbon. 




FIG. 7. (Color online) Probability densities corresponding to 
two states of the edge bands closer to Ef , near a high sym- 
metry point, for several flat ribbons. The spin polarization is 
plotted in white for spin up and blue (gray) for spin down, 
(a) Zigzag 10(1,0) ribbon, k = 0.95 TX; (b) chiral 4(7,1) rib- 
bon, k = 0.10 rX; (c) chiral 1(7,2) ribbon, k = 0.95 TX; (d) 
chiral 4(5,2) ribbon, k = 0.10 TX; (e) chiral 2(3,2) ribbon, 
A; = 0.95rX. 



In curved ribbons, edge states keep its localized char- 
acter, even for the maximum curvature. Besides the spin 
splitting (Fig. |6|, there are two main differences: the 
spatial localization length is larger for the curved ribbon 
than for the flat case, and the spin direction changes 
from the direction normal to the ribbon surface, ac- 
quiring a component in the ribbon plane, as previously 
reported.EHlHI 



VI. DISCUSSION AND CONCLUSIONS 

Our results show that the relation between the QSH 
edge states of graphene nanoribbons and both, the crys- 
tallographic orientation of the edges and the curvature 



of the GNR, allows to control the spin and spatial lo- 
calization of the ribbon edge states. Consequently, this 
may allow to achieve an efficient electrical control of spin- 
currents and spin-densities in GNRs. Taking into account 
both spin and valley degeneracy, Bloch states in graphene 
are four- fold degenerate. SOC splits them into two 
Kramers doublets and, as Kane and Mele predicted, ^^^ 
this turns graphene into a topologically non-trivial ma- 
terial. In zigzag-edge ribbons, it has already been shown 
that at each edge spin up and spin down electrons move in 
opposite directions. ^^ ^^ Since backscattering in a given 
edge requires the reversal of spin it cannot be induced by 
spin-independent scatterers. Accordingly, edge states in 
zigzag GNRs are topologically protected and hence the 
conductance of the edge states is quantized. However, 
in flat chiral-edge GNRs of finite width and for realistic 
values of the SOC strength, there is a non-zero probabil- 
ity of having electrons moving in opposite directions with 
the same spin polarization at a given edge (see fig. ItI). 
Therefore, intra-edge backscattering may occur, which 
affects the quantization of the conductance. In chiral- 
edge GNRs, spin reversal can be induced even by non- 
magnetic disorder and thus edge states do not present 
a robust conductance quantization. The appearance of 
backscattering does depend on the chirality angle, in- 
creasing for angles approaching the armchair limit. On 
the other hand, curvature breaks the inversion symme- 
try of the ribbons and Bloch states are spin-split. Elec- 
trons with the same spin and opposite propagating di- 
rections in a given edge have different energy. As a re- 
sult, backscattering is not allowed and in curved chiral- 
edge GNRs, edge states behave as robust quantum chan- 
nels. Therefore, these chiral ribbons present a magneto- 
mechanical effect: upon curving the ribbon, the spin 
channels are split in energy, thus allowing for spatially 
separated spin currents. 

Despite the weakness of SOC in these carbon systems, - 
in the scale of a few meV-, their effects in curved graphene 
and nanotubes are not negligible, due to the coupling 
of TT and a bands in curved geometries. Thus, although 
small, the effects discussed in this work may be physically 
relevant, and since in graphene the position of the Fermi 
level can be adjusted with external gate voltages, the 
control of spin currents in GNR-derived devices could be 
possible. 

In summary, we have shown that, in the presence of 
spin-orbit interaction, curvature breaks spin degeneracy 
in graphene ribbons. Flat nanoribbons with symmet- 
ric edges, either chiral or achiral, have spin-degenerate 
bands. This is due to the existence of spatial inversion 
symmetry in flat ribbons, which is broken in the curved 
cases. Furthermore, spin-orbit splitting is enhanced in 
curved ribbons due to the hybridization of the bands, 
absent in flat samples. Other mechanisms to break in- 
version symmetry, such as edge modification, are much 
less efficient to remove spin degeneracy. 

We have also explored finite-size effects in GNRs. We 
find that narrow chiral ribbons present a sizeable gap. 



despite their having a zigzag edge component, whereas 
in pure zigzag GNRs of similar width the gap is negligi- 
ble. We relate this behavior to the boundary conditions 
in chiral edges, which mix the two sublattices at each 
edge. Finally, we have studied the chirality dependence 
on the spatial localization of edge states. In narrow chi- 
ral ribbons, edge states have a nonzero density at both 
edges simultaneously, due to edge coupling. Due to the 
sublattice mixing produced by the chiral boundary con- 
ditions, edge states have a larger penetration than those 
of achiral ribbons. For wider curved ribbons, they behave 
as spin-filtered states, being localized at one edge. 
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